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The Busing-Levy IBM 704 least squares program has been applied to 
three-dimensional X-ray diffraction data from crystals of guanidinium 
aluminum sulfate hexahydrate taken with the Bond-Benedict single-crystal 
automatic diffractometer. Indications of interactions between parameters 
were evident in the early stages of refinement and were not removed in the 
subsequent cycles. Strong interactions were subsequently corroborated by 
large values of many of the correlation coefficients of pairs of parameters. 
In this case these interactions prevent refinement. The correctness of the 
general features of the structure as given in a previous paper on the gallium 
isomorph is nevertheless corroborated by the present investigation. 

To enable those who have had similar difficulties to compare results, a 
fairly detailed account is given of the course of the attempt to refine the 
structure. The effects of highly correlated parameters are emphasized. 

I. INTRODUCTION 

The purposes of the investigation to be described were manifold. An 
approximate structure of the isomorphous gallium compound has al- 
ready been reported. 1 The gallium compound with the heaviest metal 
atom among the isomorphs appeared to be best for establishing the 
general features of the structure." However, in the hope of rinding a 
closer relation between the structure and its electrical properties, it 
appeared that a refinement of the structure would be very worthwhile. 
In such a case, one would wish to have all of the atoms of more nearly 
the same scattering power; thus the guanidinium aluminum sulfate 
hexahydrate (CJ.A.S.H.) compound seemed most suitable for this pur- 

t The contribution of H. Katz to this work was made during a period of employ- 
ment at Bell Telephone Laboratories in the summer of 1959. 
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pose. Furthermore, this crystal would have the lowest linear absorption 
coefficient for all practical radiations; the importance of this feature will 
be discussed later. But probably most important, it was anticipated 
that the aluminum compound would be the one on which most measure- 
ments of various sorts would be made. This has indeed been the case. 

While our earlier paper 1 was in press, a note" 1 appeared in Kristallo- 
grafiia which gave an approximate structure for G.A.S.H. and its iso- 
morphs which differed from that reported by us. A check with our data 
indicated that the structure reported by Varfolomeeva et al* was incor- 
rect, 2 but this did not mean that the structure reported by us was neces- 
sarily correct. We had to face the question as to whether the correct 
structure might lie between the two structures or as mentioned in our 
first paper, perhaps some subtle disorder existed in the structure. In any 
case the appearance of the other result gave additional impetus to 
completion of work that had been started several years ago. 

There is a further importance of this work. The quantitative X-ray 
data were taken with the Bond-Benedict single-crystal automatic dif- 
fractometer. 4 It is the only crystal so far studied with this equipment 
and perhaps is the first X-ray structure analysis to be based on three- 
dimensional data collected automatically. Thus at least a small part of 
this paper will be devoted to an assessment of this equipment and sug- 
gestions as to future plans. 

Perhaps the most frustrating experience encountered is to find inde- 
terminate a problem which has taken considerable expenditure of time 
and effort of various sorts. One such reported problem in the field of 
X-ray crystallography is that of the determination of the structure of 
tetragonal BaTi0 3 ; this problem was found by Evans 5 to be indeter- 
minate by X-ray analysis, at the very least on the basis of the data 
collected. The results of the work on the three-dimensional data of 
G.A.S.H. indicate that the structure as originally reported by us is 
essentially correct. But we find that although a low discrepancy factor 
and standard error of fit are obtained by the least squares method of 
refinement, the structure cannot be refined; that is, convergence is not 
attained: there are parameter oscillations in each least squares itera- 
tion; some improbable interatomic distances and large error estimates 
are obtained. The cause appears to be strong interdependence of many 
of the parameters. 

In this investigation the correlation matrix is used to demonstrate the 
existence of the strongly interacting parameters. The importance of this 
approach has also been demonstrated by a recent investigation de- 
scribed in a paper written by one of us (S.G.). 
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Table I — Lattice Constants of Guanidinium Aluminum 
Sulfate Hexahydrate 



Investigators 



Wood 

Ezhkova, et al 
This work 



,i,.\ 



11.77 ± 0.04 
11.737 ± 0.002 
11.75 ± 0.02 



e,A 



8.98 ± 0.03 
8.948 ± 0.002 
8.94 ± 0.01 



II. CRYSTAL DATA 

Guanidinium aluminum sulfate hexahydrate, C(NH 2 )3A1(S04)2-6H 2 0, 
is isostructural with the previously reported 1 gallium compound. The 
morphology and unit cell dimensions have been reported by Wood. 7 
Lattice constants have also been reported by Ezhkova et al. 8 The central 
values of our lattice constants, obtained from careful measurement of 
Buerger precession camera photographs, differ from those reported in 
both of the aforementioned papers, but are in better agreementf with 
those of Ezhkova el al. For purposes of comparison, the variously re- 
ported values are listed in Table I. 

As described earlier, 1 the most probable space group to which the 
crystal belongs is P31m and the unit cell contains three formula units. 
The molecular weight of the Al compound is 387.29, the volume of the 
unit cell is 1,069 A , and the X-ray density is 1.804 g/cc. 

III. DETERMINATION OF THE STRUCTURE 

The determination of the structure has been described in the paper 
on the gallium compound. The evidence for the correctness of the general 
features of the structure described in that paper, including the orienta- 
tion of the guanidinium ions, is conclusive as will be shown subsequently. 



IV. EXPERIMENTAL 

The Bond-Benedict single-crystal automatic diffractometer 4 was used 
to collect the three-dimensional data. Some changes from the original 
design of the instrument and in the electronics were made before the 
final data were taken. A detailed description of these changes must be 
left, to the original authors. However, it should be mentioned that for 
these particular data (which were taken in 1956), a proportional counter 
replaced the Geiger counter and the "back-set" correction 4 was virtually 

t Dr. E. A. Wood and Mrs. V. B. Compton have informed us that their recent 
measurements of lattice constants of (i.A.S.H. give values which agree more closely 
with those of Ezhkova et al. H and of the present work. 
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eliminated by circuitry changes. Also, the internal geometry of the 
collimator was changed to square cross section. 

The need for a collimator with square cross section derived from the 
mechanics of the instrument. The "back-setter" produces a jarring of 
the goniometer head which could at times translate the crystal very 
slightly out of the original alignment in the X-ray beam. If the beam 
has a circular cross section, slight deviation from coincidence of crystal 
cylinder and rotation axes causes significant differences in intensity 
when the diameter of the crystal is large relative to the beam cross 
section. This is not true of a beam with a more or less square cross 
section. 

Of course, one would not have to worry too much about this if small 
crystals were being used. However, for this instrument and the use of 
the usual type of sealed X-ray tube, it is necessary to use large crystals 
to obtain the data. (This will be discussed further later.) 

Two cylindrical crystals were used to obtain the data attainable by 
this instrument with CuKa radiation and a pentacrythritol mono- 
chromator. The crystal aligned along the c-axis had a diameter of 0.67 
mm; the crystal aligned along the [20- 1J direction (orthohexagonal 
A -axis) had a diameter of 0.54 mm. With a linear absorption coefficient 
for CuKa radiation of 48.7 cm -1 , the values of nR for these crystals are 
1.64 and 1.32 respectively. 

As described in the paper by Bond, 4 the single-crystal automatic 
diffractometer works on a principle similar to that of the equi-inclination 
Weissenberg camera. With CuKa radiation, seven levels were obtain- 
able about the c-axis and fifteen about the orthohexagonal A -axis. 

Data from a particular level n were collected as follows: The align- 
ment of the crystal was checked. This was done in two ways whenever 
possible. A microscope could be used to align the crystal cylinder axis 
with the rotation axis of the instrument. The equi-inclination angle was 
calculated and the crystal set to this angle. The arrangement of the 
counter of the instrument is always set so that the diffracted beam is 
incident perpendicularly to the window. Thus the counter is actually 
moved to twice the angle of the crystal from the zero level situation. 
If a particular reflection (for example, 00-/ on the Ith level about the 
c-axis) was observable when the counter angle was equal to zero de- 
grees for a given layer, this reflection was used to readjust crystal and 
counter. 

To obtain the weak intensities, the diffraction unit settings were usu- 
ally 40 kv and 20 ma. To obtain the stronger reflections, proper settings 
of the voltage and tube current were made so as to record enough 
moderate reflections to establish a scale between the two patterns. 
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Integrated intensities, crystal angles and counter angles for each level 
were recorded automatically by the Leeds-Northrup two-pen recorder 
as described in the papers by Bond and Benedict. 4 As indicated above, 
resetting was made manually for each new level. 

Following the collection of the data by the recorder, it was necessary 
to index the data: This was the most time-consuming (i.e., on a man- 
hour basis) part of the data processing required to obtain the observed 
amplitudes. The indexing was carried out with the use of the plotting 
device. (The indexing problem will be discussed further later.) 

Following the indexing of all the data, the usual absorption, Lorentz- 
polarization and Tunell 9 rotation factorsf were applied to extract the 
relative | F |\ (The polarization correction is for monochromatizcd 
radiation.) The calculation was programmed for the IBM 704 by R. G. 
Treuting. The corrections calculated were based on the formulae J 
given by Bond and the tables used for the absorption corrections are 
those given in Bond's paper. 10 The program written by Treuting put the 
resultant | F„ | 2 's or | F„ |'s out on cards as well as on a print-out. The 
individual Lorentz-polarization, absorption and Tunell rotation factors 
were also printed out for each reflection for each layer on which it 
appeared. 

Having extracted the | F„ | 2 's for each layer about each of the two 
axes, the next step involved an iterative cross-calibration process to 
bring the values to the same basis. An IBM 704 program written by 
W. R. Romanow allowed us to apply constant factors to the sets of 
| F„ | 2 put out by the intensity correction program. Romanow's program 
also put out new cards so that we could apply a different constant to the 
new values if necessary. 

When we felt we had arrived at the best values, it was decided to 
carry out the subsequent least squares refinement on the basis of the 
| F | values. Using a short program written by Romanow, square roots 
were taken of all the | F | 2 's and put out on cards. Those that came from 
layers about the orthohexagonal -4 -axis were then sorted on the values 
of / for ease in setting up the data for the least squares refinement. 

As described in the Bond-Benedict papers, some reflections do not 
get entirely into the counter; thus, in order to be sure that all are ob- 



t The proportional counter employed had a linear response to counting rates of 
over 20,000 cps. Because for even the strongest reflections, observed counting rates 
over 10,000 cps gave integrated intensities which went off scale on the recorder, no 
dead-time correction 4 was necessary for any of the reflections. 

J The formula for Pi on p. 380 of Bond's paper should read 

P L = T sin 20 / 1 + ( ? ~ *'"' -- ) (i + C os 16V - j-v~ (l + cos 29) 

I \ \0 + $ eos v l + g 
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tained, the instrument was designed to obtain each reflection twice. 
For this reason the counter has a 4° window. Even at that, not all the 
reflections of a given form will have the same intensity, but usually 
about a twofold axis, a form of reflections of moderate intensity will 
have two with the same intensity. About a threefold axis, perhaps eight 
of twelve reflections from a given hk • I form will have the same intensity 
or 12 out of 16 of a given hk-Q form. Unfortunately, the weaker reflec- 
tions do not give as good results as the moderate to strong ones. In the 
case of the c-axis layers, if there was a variation in the height of peaks 
which appeared to have been fully in the window, the value taken for 
the integrated intensity was the average of the several peaks. In the 
case of the orthohexagonal A -axis layers, because there were fewer 
peaks contributing to a form and therefore a greater possibility that 
only one peak was squarely in the window, the value recorded in most 
cases was the measure of the highest peak. 

In taking the averages of observed structure amplitudes, the weighting 
was in accord with the above. For example for a given | F hk .i |, h,k,l ?± 0, 
the value from the c-axis layer was weighted four times and a value from 
an orthohexagonal A -axis layer once. The standard deviation was cal- 
culated in accordance with the analysis given in Chapter 16 of the book 
by Dixon and Massey 11 and as suggested earlier by Ibers. 12 However, 
for the unobserved, the standard deviation was taken as equal to half 
the minimum observable. For | F 00 .i |'s which would have unity weight 
since they appear only once, the a was taken in accordance with a 
subjective estimate comparing the particular | F 00 .i | with others of 
similar value. The agreement between or among | F n |'s from the same 
form but from different layers was quite good generally except for the 
weakest reflections. 

In the CuKa sphere, there is a total of 895 X-ray forms of guanidi- 
nium aluminum sulfate hexahydrate. The geometry of the Bond diffrac- 
tometer allows us to observe only 842 of these. Of those possibly ob- 
servable by the instrument, only 546 were actually observed. 

V. ATTEMPT TO REFINE THE STRUCTURE 

Because the major point of this paper is to demonstrate that the 
refined structure under discussion is effectively unattainable from the 
X-ray diffraction data, it seems worthwhile to give some of the details 
of the calculations. To make such a discussion simpler, the pertinent 
data are collected in tables. In Tables II and IV the values of parameters 
and some other important information are listed. In Table II two col- 
umns are assigned to each cycle; the left one lists the starting parameters, 



Table II — Results of Least Squares Calculations (Fihst Set op Weights) 



Cycle 4 



Scale 1 


0.6667 


0.553 




0.583 


0.550 


0.542 




0.5395 






0.5612 




0.5552 


0.5467 


0.5552 


0.5569 


0.5557 










Fictori 2 


0.6667 


0.51)2 
0.51)8 




0.569 


0.540 


0.534 




0.5552 






0.5889 




0.6162 


0.6200 


0.6162 


0.6161 


0.6186 










I 


0.6667 




0.566 


0.550 


0.581 




0.5695 






0.5829 




0.5879 


0.5703 


0.5679 


0.5890 


0.5933 










0.6667 


0.580 




0.619 


0.580 


0.580 




0.5879 






0.6221 




0.6327 


0.6132 


0.6327 


0.6208 


0.6161 










5 


0.6667 


0.579 




0.564 


0.580 


0.619 




0.6187 






0.6411 




0.6492 


0.6909 


0.6492 


0.6396 


0.6449 










6 


0.6667 


0.610 




0.575 


0.610 


0.630 




0.6494 






0.6782 




0.6887 


0.6978 


0.6887 


0.7116 


0.7327 










7 


0.6667 


0.567 




0.5B0 
0.1)93 


0.570 
0.430 


0.609 




0.6533 






0.6960 




0.6887 


0.6739 


0.6887 


0.6881 


0.6854 










B 


0.6667 


0.1)33 




0.545 




0.6096 






0.6740 




0.6826 


0.6873 


0.6826 


.0.6836 


0.6856 










Atom 1 x 


0.553 


0.551 




0.554 


0.551 
0.407 


0.558 




0.5624 






0.5632 


0.5540 


0.5600 0.5333 0.5626 


-0.5533 


-0.5531 

-0.3349 


-0.5522 




-0 


.4630 
.2740 


-0.4614 


H (ll) 1 


0.418 


0.1(07 




0.390 


0.409 


0.3333 


0.3379 






0.3344 


0.3333 


0.3386 0.3333 0.3312 


-0.3333 


-OJ338 







0.2791 


i 


0.520 


0.527 


, 


0.540 


0.527 


0.532 




0.5544 






0.5567 




0.5492 


0.5639 


0.5492 


0.5350 


0.5275 


0.5236 





.535 


0.5321 


B or P,, 


2.00 


6.13 


4.50 


10.2 


0.01086 


0.01611 


0.01000 


0.01173 






0.00887 














0.02497 


0. , 


.nrtRQn 




5k 










0.01086 


0.02353 


0.01000 


0.01478 






0.00040 


0.00700 












0.01493 


0- 00700 




K33 










0.01408 


0.02927 
0.00412 


0.02000 


0.01223 






0.01817 














0.02670 





.01820 




gl2 










0.00543 




0.01032 






0.00503 














0.01499 





.00500 




B 23 










0. OOOOO 
0.00000 


0.00650 
0.00135 




0.00188 

0.00084 






O.0O516 
0.00506 














0.01027 
0.00615 






.00520 

.00410 




Atoa 2 x 


0.31)2 


0.339 




0.345 


0.339 


0.340 




0.3422 






0.3435 




0.3437 


0.3459 


0.3437 


0.3451 


0.3449 








0.3455 


o (v) y 


0.120 


0.113 




0.116 


0.113 


0.114 




0.1166 






0.1252 


0.1200 


0.1191 


0.1176 


0.1191 


0.1165 


0.1166 








0. 1159 


z 


-0.220 


-0.221 




-0.217 


-0.221 


-0.222 




-0.2114 






-0.2010 




-0.2064 


-0.2075 


-0.2064 


-0.2008 


-0.1972 


-0.1976 
0.00843 






-0.2009 


B or P., 


2.00 


1.75 


1.50 


0.74 


0.00362 


0.00614 




0.00748 






0.01087 


0.00850 


















222 

s 










0.00362 


0.00314 




0.00219 






0.00215 














0.00307 
















0.00469 


0.01295 
0.00563 




0.01947 






0.01511 














0.02193 
















0.00181 


0.00400 


0.00510 


0.0*350 


0.00611 


o.oo4oo 












0.00375 
















0.00000 


0.00221 




-0.00016 






-0.00338 














-0.00333 
















0.00000 


0.00029 




-0.00022 






-0.00222 














-0.00206 








Aid 3 x 
O(VI) y 

B or !,, 


-0.314 

-0. 120 

0.258 

2.00 


-0.321 

-0.120 
0.2»0 
1.90 


1.75 


-0.333 
-0.126 

0.236 
1.35 


-0.321 
-0.120 

0.240 

0.00423 


-0.324 
-0.120 

0.235 

0.00271 




-0.3222 
-0.1171 

0.2349 
0.00080 






-0.3235 
-0. 1 130 

0.2423 
-0.00106 


-0.1170 
0.00100 


-0.3214 -0.3208 
-0.1122 -0.1.150 -0.1119 
0.2362 0.2400 0.2328 


-0.3214 
-0.1150 
0.2400 


-0.3212 

-0.1139 
0.2448 


-0.3211 
-0.1137 
0.2450 


0.00287 

0.00149 
0.01905 
0.00066 
-O.OOOBl 
-0.00057 






-0.3218 
-0.1173 

0.2439 


?33 










0.00423 
0.00547 


0.00033 
0.01067 


0.00100 


0.00094 
0.01162 






0.00344 
0.01523 




















2l2 










0.00211 


O.OOI53 




-0.00078 






-0.00155 




















2l3 










0.00000 


-0.00002 


0.00000 


0.00028 






-0.00181 




















»23 










0.00000 


-0.00030 


0.00000 


-0.00052 






0.00012 




















Mom 14 x 

o(ix) y 

B or | (| 


-O.II77 
0.333 

-0.066 
2.00 


-0.1)71 
0.330 

-0.058 
0.73 


1.50 


-0.475 
0.332 
-0.059 
-0.22 


-0.471 
0.330 
-0.058 

0.00362 


-0.471 

0.327 

-0.059 

-0.00259 


0.00000 


-0.4687 
0.3258 
-0.0593 
-0.00104 





.00100 


-0.4673 
0.3263 

-0.0577 
0.00063 




-0.4654 
0.3288 
-0.0588 


-0.4702 
0.3261 
-0.0654 


-0,4654 
0.32BB 
-0.0588 


-0.4649 
0.3270 
-0.0599 


-0.4654 
0.3272 
-0.0546. 


-0.0520 






-0.4677 
0.3291 
-0.0504 


l?2 

§33 










0.00362 
0.00469 


0.00186 
0.01372 




0.00003 
0.01363 






0.00036 
0.01626 














0.00266 
0.00009 
0.01613 
0.00127 
-0.00122 
0.00174 





.00120 




Si? 










0.00181 


1-0.00006 


0.00000 


-0.00223 





.00000 


0.00017 




















bi 










0.00000 


-0.00052 


0.00000 


-0.00083 






0.00083 




















"h 










0.00000 


-0.00175 


0.00000 


-0.00028 





.00000 


0.00147 




















Atom 5 x 
(X) y 

B or p,, 


O.u.76 
-0.336 
0. 174 
2.00 


0.1)65 
-0.31)3 
0.172 
2.1)7 


2.00 


0.459 

-0.347 

0.176 

2.16 


0.465 
-0.343 
0.172 
0.00483 


0.466 
-0.340 
0.174 
0.00924 




0.4614 
-0.3391 
0.1775, 

0.00444 






0.4575 
-0.3380 
0.1816 
0.00459 




0.4635 
-0.3417 
0.1792 


0.4636 
-0.3389 
0.1762 


0.4635 
-0.3417 
0.1792 


0.4664 
-0.33BB 
0.1717 


0.4647 

-0.3405 
0.1756 


-0.3391 
0.1758 






0.4665 

-0.3398 
0.1810 


g22 










0.00483 


0.00118 




0.00097 






0.00385 














0.00099 
0.00142 





.00200 




£ 3 3 










0.00626 
0.00261 


0.00588 




O. 00905 






0.01667 




















p 13 










0.00112 




-C. 00097 






0.00037 














0.01636 
















0.00000 


-0.00422 


0.00420 


-0.00243 






-0.00322 














0.00135 








p» 










0.00000 


0.00060 




-0. 00062 






-0.00187 














0.00007 
-0.00150 








Atom 6 X 
N (1) z 
B or B , 


0.113 
0.500 
2.00 


0.113 

0.1)97 

-2.50 


3.00 


0.106 
0.519 
-4.60 


0.113 

0.497 , 
0.00724 


0.112 

0.500 

-0.01095 


0.00500 


0.1213 
0.5037 

0.00378 






0. 1258 
0.5068 
0.01679 


0.1130 

0.0038 
0.01480 


-0.1130 -0.1 138 
0.5038 0.4930 


0.1130 
0.5038 


0.1135 
0.5177 


0.1132 
0.5065 


0.5175 






0.1198 
0.5137 


gas 










0.00938 


-0.00125 


0.00500 


0.01479 






0.00380 












0.00798 





.00380 
.01480 




lis 










0.00362 


-0.00812 


0.00000 


0.00928 





.0032 


0.02214 


0.00320 












-0.00095 
0.01743 
-0.01087 







"is 










0.00000 


-o.ooe39 


0.00000 


-0.00028 






-0.00593 


0.00000 












o.oaiao 
0.00000 




Alom 7 X 
(VIII) z 
B or P,, 


-CM 
-0.135 
2.00 


-0.135 

-0.121) 

1.53 


1.50 


-0.133 

-0.118 
-0.04 


-0.135 
-0.124 
0.00362 


-0.134 
-0. 125 

0.00529 




-0. 1360 
-0. 1208 
0.00732 






-0.1387 
-0.1149 

0.00448 
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the right, the calculated "corrected" parameters. A blank space in the 
left column indicates that the last previous calculated value was the 
starting value for the particular parameter. In the cases of cycles 9 and 
10, all of the parameters had the last previous calculated values of 
cycles 8 and 9 respectively. 

The order in which the atoms are listed in Tables II and IV is not 
the same as that of the paper 1 on the gallium isomorph, but the atom 
labeling is. In writing the special position symmetry patch for the 
Busing-Levy 13 IBM 704 least squares refinement program, it is most 
convenient to list the atoms in general positions first. Then to avoid 
mistakes in the listing of results, it is best to leave the order the same 
as that of the output of the program. 

In the calculation of structure amplitudes the following atomic scat- 
tering factors were used: for O, Al 3+ , X and C, those of Berghuis ec al; 1 * 
and for S, those of Viervoll and Ogrim. 10 

In cycles 1 and 2, 895 reflections, all those representing independent 
forms and observable in the CuKa sphere, but perhaps not observable 
with the single-crystal diffractometer, were included. Eight of the 
parameters were scale factors, all of which were initially equal to 0.(>fio7, 
one for each value of / from to 6 and the eighth value for all the re- 
maining / values. Also in the first two cycles, isotropic temperature 
factors were used despite the fact that it was obvious that the thermal 
motions of the atoms in this crystal must be highly anisotropic. 

The starting structural parameters for the first cycle were those given 
for the gallium isomorph 1 except for changes in the S and Al tempera- 
ture factors and the //-parameter of X(II), which was inadvertently 
taken as 0.418 instead of 0.333. Xow it may be seen in Table II under 
cycle 1, that this //-parameter did not change as radically as one might 
have hoped, in fact as one might have expected, for an incorrect parame- 
ter. But the temperature factor of the atom did increase considerably, 
perhaps indicating that the atoms did not want to be at the positions 
indicated. On the other hand, the temperature factor of the N in the 
special position decreased considerably to a negative value as if to 
compensate for the other. This, in retrospect, was already indicative of 
strong interaction between the thermal parameters of these two atoms. 
Another important change was the large one, to —0.392, in the value of 
the O(III) z-parameter; this implies a very short S — O distance, 1.31 
A, in one set of the S() 4 groups. 

The estimated error of fit '''•'' at the end of the least squares calculation 
of cycle 1 was very much lower than the first computed error of fit, 13 
and it appeared that by readjustment of some of the temperature factors 
we could go a step further toward convergence before changing to aniso- 
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tropic temperature factors. Initially cycle 2 showed that even with the 
readjustment of temperature factors, the R value f had dropped from 
0.473 to 0.303, the weighted R from 0.299 to 0.193. But the error of fit 
was higher than that estimated in cycle 1 on the basis, of course, of the 
parameters computed in that cycle, some of which were physically 
impossible (i.e., negative temperature factors). 

However, cycle 2 ended with an estimated error of fit somewhat lower 
than that of cycle 1. The N(II) //-parameter decreased toward the value 
which we believe to be the more nearly correct one, but the N(II) B 
value increased greatly and the N(I) B value became a large negative 
value. Also the ^-parameter of X(I) decreased to imply an unlikely short 
C— N distance. Changes in the S and Al positional parameters were not 
laige but seveial oscillations occurred. The O(III) (atom 10) ^-parameter 
letuined to -0.400, but even this value implied a rather short S— O 
distance, 1.37 A. 

At this point, it seemed necessary to change to anisotropic thermal 
parameters. The Busing-Levy program will compute these from the 
isotropic thermal parameters using the following relations: p n = Ba* 2 /A; 
fi l2 = (#a*b*cos7*)/4;etc. 

The starting parameters were those computed in cycle 1 and adjusted 
for cycle 2 (see Table II). For cycle 3, a critical estimate of the reflec- 
tions really observable by the single-crystal automatic diffractometer 
was made. This resulted in the removal from the calculation of 43 
unobserved reflections, some of which had rather high calculated struc- 
ture amplitudes when compared with the respective estimated threshold 
values. Included in cycle 3 was a rejection test: that is, when A/ a was 
> 10.00, the reflection was not counted in the calculation of the R 
values or the standard error of fit, nor was it included in the least squares 
calculation. This reduced the number of *Wa used in the least squares 
calculation to 790. (Unfortunately the R values and the calculated 
amplitudes computed in this cycle have been lost. ) 

The estimated error of fit resulting from the cycle 3 least squares 
calculation decreased from 4.99 to 2.30, an apparently tremendous im- 
provement. However, the still incorrect N(II) //-parameter did not 
improve; also the values of the N(H) thermal parameters greatly in- 
creased. The O(III) values still implied a short S— O distance. The 
C(I) z-paramcter indicated possible nonplanarity of the guanidinium 
ion in the special position, but this parameter also had an apparently 



t Unless otherwise stated, the K value is that for the independent F hk .i's, i.e., 
multiplicity is neglected. This is the R value calculated by the Busing-Levy pro- 



gram. 
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large estimated error, 0.0115, indicative of potential difficulty. Twelve of 
the atoms had calculated thermal parameters which did not satisfy all 
the criteria for physical reality (see Ref. 13). Therefore, for cycle 4 some 
of the thermal parameters had to be adjusted to satisfy these criteria. 
Also, the N(II) //-parameter was corrected. The R value and error of fit 
decreased considerably since cycle 2, but the weighted R value increased 
slightly. The same rejection test as used for cycle 3 allowed 809 reflections 
to be included in the cycle 4 calculation. The least squares calculation 
led to an estimated error of fit of 2.23, not too different from that esti- 
mated in the previous cycle. 

In cycle 4, the values of the N(II) thermal parameters decreased, 
indicating that the high values had been caused by the wrong //-parame- 
ter. One would prefer to think, however, that the //-parameter should 
have tended to approach the correct value rather than to have the 
thermal parameters act as if the atom should be removed. This time the 
^-parameter of N(I) (atom 6) became rather large, implying too large a 
C — N distance. A number of the other positional parameters showed 
oscillation, and again there were twelve atoms which had thermal 
parameters not satisfying the criteria for physical reality (Table II). 
The O(III) ^-parameter continued to imply a short S — distance. 
The C(II) and N(II) atoms did not have the same values in 2-parame- 
ter, nor did the C(I) and N(I) atoms have the same 2-parameter. Also, 
in this cycle many of the scale factors, especially s 8 , had almost reached 
their starting values after having at first decreased substantially. 

The necessary adjustments were made on the thermal parameters 
before cycle 5 was carried out. Also, the rejection test was removed. Five 
reflections which appeared to have substantial contribution from the 54 
hydrogen atoms or to have suffered from extinction were given zero 
weight. Thus, of the 852 reflections, 847 were used in the cycle 5 least 
squares calculation. Because some of the initially estimated o-(F„)'s were 
really very small, a few of these also were readjusted. Initially the R 
value was 0.198, while the weighted R decreased to 0.139, this latter 
reduction resulting mostly perhaps from the few adjustments made on 
the o(F„)'s. The error of fit for the 847 reflections was larger than for 
the 809 of the previous cycle. The calculated estimated error of fit after 
the least squares calculation did decrease, however. 

But in cycle 5 there was no improvement in the way the calculation 
was going. There were further oscillations, and, very important, the 
0(H) — X(II) distance continuing from cycle 3 was short, whereas the 
C(I) — N(I) distance continued to be long. Considering the guanidinium 
ion to be planar, the C — X distances were respectively 1.22 and 1.48 A, 
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the average is 1.35 A in good agreement with the acceptable guanidinium 
C — N value 1.34 A. 16 Again this indicated interaction between the N(II) 
x- and (/-parameters and the N(I) x-parameter. Also, the parameter 
values of the S(I) and O(III) atoms still indicated an improbably short 
S — O distance. There were other indications of interaction: for example, 
the (/-parameters of the O(V) and O(VI) atoms (2 and 3 respectively) 
behaved strangely, that of O(V) indicating an improbably large [S0 4 ] 
0—0 distance, that of O(VI), too small an [S0 4 ] 0—0 distance. 

It seemed at the time, however, that there might be other possibilities 
for explaining the course of events in the attempt to refine the structure. 
For example, there could be many reflections to which the hydrogen 
atoms would contribute, and, perhaps particularly because this is a non- 
centrosymmetric structure, the affected structure amplitudes were hav- 
ing a detrimental effect. Therefore, in cycle 6 all reflections for which 
sin 2 /X 2 < 0.0800 were given zero weight. Necessary adjustments were 
made in thermal parameters (Table II); the N(I) and N(II) positional 
parameters were readjusted each to yield the C — N distance 1.34 A; 
and the 0(V) and O(VI) (/-parameters were adjusted to yield more 
reasonable [S0 4 ] 0—0 distances. The R value for the 755 amplitudes 
(with nonzero weights) was 0.200, weighted R = 0.128 and error of fit, 
2.82. 

In the cycle 6 least squares calculation, only 43 parameters were 
varied: the scale factors and all positional parameters except the N(I) 
re-parameter. The estimated error of fit decreased to 2.38, but this cycle 
was also discouraging in that again there were oscillations and some 
rather large changes in parameter. The S(I) — O(III) distance continued 
to remain improbably short; the O(VI) (/-parameter again implied too 
short an [S0 4 ] O—O distance; and the values of the N(II) x- and y- 
parameters implied a C(II)— N(II) distance of 1.25 A. 

In the paper on the gallium isomorph, 1 we had concluded that the 
arrangement of the guanidinium ions on the threefold axes were related 
to that at 3m to close approximation by J,|,0 and f,£,0 - (u,0,w; 
0,u,w; u,u,w). However, some doubt remained, and therefore it was 
decided to try some different orientations of the guanidinium groups. 

For cycle 7, the N(II) parameters were readjusted, presumably back 
to the starting parameters of cycle 6. However, a card-punch error 
(0.5333 instead of 0.5533) was made in the s-parameter. The N(I) 
parameter was set to -0.1130. This we shall call the ( — ,— )t orienta- 

t This symbolism is derived as follows: The ± orientations of N(I) are those for 
which in (x,0,z) of positions 3r, zmo) = ±M where u is very nearly +0.113. The ± 
orientations of N(II) are those for which in (x,y,z) of positions Gd, Xkiu) = i ± », 
y = §. Thus (-,-) here means that x n ,d = - 0.113, jniid = 0.220, t/niii) = 0.667. 
By symmetry the latter two are equivalent to 0.553 and 0.333 respectively. 
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tion. The positional parameters of O(VI) were also readjusted. The R 
value for the 755 reflections increased to 0.250, the weighted R to 0.187, 
and the error of fit to 4.10. In cycle 7 all scale and positional parameters 
were varied. At the end of the cycle, the estimated error of fit was 3.53. 
The C(II) — N(II) distance again was too short, ~1.21 A; again the 
()(VI) //-parameter decreased from the adjusted value; the difference in 
the C(I) and N(I) z-parameters increased. Also again there were oscil- 
lations. The results of cycle 7 did not look promising. 

In cycle 8, the (+,+) arrangement of the guanidinium ions was tried 
with the other starting parameters the same as those used in cycle 7. 
In this case the R value for the 755 amplitudes was 0.231, weighted R, 
0.155, and error of fit, 3.40. Again only scale and positional parameters 
were varied. The estimated error of fit obtained at the end of the least 
squares calculation was 3.14. The results of this cycle looked promising. 
The C — N distances looked good; the O(V) and O(VI) parameters 
were not too bad. However, the S(I) — O(III) distance still looked 
improbably short, The agreement for individual amplitudes actually did 
not look as good as it did in cycle 6, but it was felt that perhaps some of 
this poorer agreement resulted from hydrogen contributions and/or from 
required changes in thermal parameters. 

It was decided to continue to cycle 9 using the values of scale and 
positional parameters obtained in cycle 8. The R value increased to 
0.240; the weighted R value decreased to 0.140; the error of fit was very 
close to that previously estimated. Despite this, the parameter results 
of this cycle (Table II) looked even better than those of the previous 
cycle, but the S(I) — O(III) distance continued to be improbably 
short. 

The scale and positional parameters resulting from cycle 9 were used 
in cycle 10. There was not much change in R, weighted R or error of fit. 
In cycle 10, all scale and positional parameters which had changed less 
than la in cycle 9 were held constant and all thermal parameters were 
allowed to vary. The estimated error of fit at the end of the cycle was 
2.55. It appeared that the thermal parameters of the N(II) atom in- 
creased considerably as if trying to eliminate this atom, and as before 
this seemd to be an indication that the N(II) atom was not placed 
correctly. Also as if to compensate, the previously large /3 33 of N(I), 
0.01480, decreased to -0.00095. Eight of the atoms had thermal 
parameter matrices which were not positive definite. 

With this continued disappointment, another notion became more 
important. Was it possible that the structure given by Varfolomeeva 
et al 3 was correct? It seemed advisable to make the calculation with the 
model proposed by those authors. The results proved that the structure 
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cannot possibly be correct, The initial R was 0.559, weighted R, 0.473 
and error of fit, 10.38 for the 755 reflections. Examination of the calcu- 
lated and observed amplitudes showed a great many very large dis- 
crepancies indicative of an improbable structure. Only the scale and 
positional parameters were varied in the least squares calculation. 
Thermal parameters for the N atoms were those initially used in cycle 6. 
All other thermal parameters were essentially those obtained in cycle 10 
with necessary adjustments made. The initial and final positional 
parameters are shown separately in Table III. The estimated error of fit 
was 8.92, indicating no real possibility of convergence. The parameter 
changes were mostly drastic. The N(I) .r-parameter, for example, would 
imply a C(I) — N(I) distance of 1.16 A. Interestingly enough, the 
S(I) — O(III) distance continued to remain very short. 

In cycle 12, the guanidinium ions on the two three-fold axes (i.e., at 
f f and f f ) were turned 30° from their original positions. The thermal 
parameters were the same as those used initially in cycle 11 and are 
shown in the next to the last columns of Table II. The R value was 
0.238, weighted R, 0.154, and error of fit, 3.38 (the latter two being 
somewhat higher than for the starting parameters of cycle 10). The 
estimated error of fit obtained from the least squares calculation was 3.14. 
The results of this calculation did not look promising. The C(I) — N(I) 
distance was large; there was an extraordinarily large change in the 
2-parameter of O(VIII). Also, agreement of many individual amplitudes 
was poorer than for the very first orientation of the guanidiniums. In 
fact, from the calculations of cycles 7-10 and cycle 12, it had become 
apparent that the ( + ,— ) orientation was indeed the best. It also ap- 
peared that disorder or rotation f of the guanidinium ions was highly 
unlikely unless very subtle. In the case of complete disorder or the 
equivalent free rotation, there would be no contributions from the nitro- 
gen atoms to the amplitudes F hk . h h — k ^ Zn, exactly as in the case 
of the ( + , + ) orientation. This alone makes it appear that the orig- 
inally reported 1 ( + ,- ) orientation of the guanidinium ions was corrob- 
orated. 

In cycle 12, the normal equations and inverse matrices were obtained. 13 
Examination of the inverse matrix sho wed that there were large values 
of correlation coefficients, Plj = btj/'s/bi&jj , for many pairs of parame- 
ters. A few examples are : 

f Two reports 1718 based on nuclear magnetic resonance investigations of G.A. 
S.H. mention the possibility of rotation of the guanidinium groups. We have 
learned (by private communication) from, and have been permitted to quote, the 
author, D. W. McCall, of one of these,' 7 that further investigation now indicates 
that this rotation is highly unlikely. 
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2o(V) — 2 0(II)j 0-81 

^o(v) — z o(i)i 0.58 

^o(ix) — 2 (vii), 0.84 

^o(iv) ~~ ^o(iri)> 0.65 

«s(i) - Zsui), 0.96. 

It is noteworthy that the correlation coefficient for x Nln) — £ N(I) was 
very low, 0.10; it will be seen later that this low value resulted from the 
incorrect orientation of the guanidinium (II) ions. 

It seemed unlikely that the weighting scheme could be the cause of 
the difficulties encountered. Nevertheless, it was decided to try a 
weighting scheme radically different from that used in the first twelve 
cycles. 

In cycle 1' (Table IV), all amplitudes with sin 2 0/X 2 < 0.0800 were 
still weighted zero. Also all unobserved amplitudes were to be weighted 
zero and all observed, unity. However, a number of amplitudes which 
should have been weighted zero were weighted unity, and a few which 
should have been weighted unity were weighted zero. This left 534 
reflections included in the least squares calculation. The initial parame- 
ters were those from cycles 9 and 10, except for the N's which were 
started at the exact ( + ,— ) orientation and the O(III) ^-parameter 
which was started at —0.405 to give an S — O value closer to 1.48 A. 
The R value was 0.204, weighted R, 0.149 and error of fit 2.19 for the 
534 amplitudes and these parameter values. The least squares calcula- 
tion gave an estimated error of fit of 1.90. Again the S(I) — O(III) distance 
decreased to 1.38 A, the C(I) — N(I) distance increased again and the 
C(II) — N(II) distance decreased again. Some of the other distances are 
listed in Table V. 

Starting with this calculation, the vector Vi — Z(\/wDi)(\/wA) was 
obtained as outputf as well as the direct and inverse matrices, 13 the 
purpose being to see whether Ap/s from the diagonal term approxima- 
tions would be much different from those obtained by the exact solution 
of the normal equations. Not many of these were checked in this and 
subsequent cycles, but enough differences were found to indicate the 
importance of the off-diagonal terms. 

It appeared that it would be most convenient to have the correlation 
or normalized inverse matrix to examine in each cycle. A program patch 
to enable us to do this was written by Misses D. C. Leagus and B. B. 
Cetlin. 



t The program patch for this calculation was written by Miss D. C. Leagus. 



Table IV — Results of Least Squares Calculations (Second Set of Weights) 







Cycle 1» 


Cycle 2' 


Cycle 3' 




fiZfl 


le 4! 






Cycle 5' 








Paramet 


nrs 












a_ 


b_ 




a_ 


b 


c 


d 


i J 


Scale 


1 


0.5557 


0.5846 


0.5627 


















0.5420 


0.0079 


Factors 


2 

3 
It 


0.6106 
0.5933 


0.6058 
0.5951 


0.6069 
0.6032 


















0.5620 
0.5472 


0.0086 
0.0072 




0.6161 


0.5841 


0.5903 


















0.5519 


. 0092 




5 


0.6449 


0.6503 


0.6539 


















0.5724 


0.0136 




6 


0.7327 


0.6826 


0.6849 


















0.6089 


0.0208 




7 


0.6854 


0.6689 


0.6760 


















0.5779 


0.0279 




8 


0.6856 


0.6693 


0.6751 


















0.5649 


0.0427 


Atom 1 


X 


0.2200 


0.2211 


0.2184 0.2237 


0.2200 0.2171 




0.2148 


0.2155 




0.2169 


0.2169 


0.2145 


0.2176 


0.0027 


N (II) 


y 


-0.3330 


-0.3377 


-0.3399 


-0.3330 -0.3386 




-0.3412 


-0.3415 




-0.3147 




-0.3419 


-0.3432 


0.0026 




z 


0.5350 


0.5465 


0.5521 


0.5579 




0.5552 


0.5559 




0.5591 
0.00424 


0.5590 
0.00422 


0.5580 


0.5581* 


0.0044 




ft 22 
a 33 
5l2 


0.00890 






0.00229 




0.00431 


0.00586 






0.00456 


0.00255 




0.00700 






0.00970 




0.01172 


0.01386 




0.01191 


0.01222 




0.01221 


0.00330 




0.01820 






0.02051 




0.02269 


0.02108 










0.01509 


0.00470 




0.00500 






0.00227 




0.00345 


0.00427 




0.00494 


0.00508 




0.00571 


0.00201 




P 23 


0.00520 






0.00293 




0.00185 


0.00170 




0.00228 


0.00226 




0.00191 


0.00233 




0.00410 






0.00039 




0.00181 


0.00217 




0.00299 


0.00309 




0.00292 


0.00269 


Atom 2 
(V) 


X 


0.3449 


0.3412 


0.3429 








0.3419 




0.3430 


0.3428 


0.3423 
0.1140 


0.3428 


0.0025 


y 


0.1166 


0.1178 


0.1167 








0.1201 




0.1174 


0.1173 


0.1165 


0.0022 




z 


-0.1977 


-0.2072 


-0.2076 








-0.2043 




-0.2005 


-0.2002 


-0.2008 


-0.2056 


0.0043 




fill 


0.00843 












0.00741 




0.00689 






0.00517 


0.00248 




S 22 

233 


0.00307 












0.00010 


0.00300 


0.00335 






0. 00244 


O.OOI73 




0.02193 












0.2388 




0.02161 






0.00999 


0.00473 




8' 2 

P 23 


0.00375 












0.00267 




0.00396 


0.00402 




0.00357 


0.00153 




-0.00333 












-0.00262 




-0.00054 


-0.00085 




-0.00060 


0.00230 




-0.00206 












-0.00217 










-0.00255 


0.00207 


Atom 3 


X 


-0.3211 


-0.3225 


-0.3243 








-0.3215 




-0.3233 
-0.1149 


-0.3231 


-0.3230 


-0.3232 


0.0022 


(VI) 


y 


-0.1137 


-0.1155 


-0.1162 








-0.1143 




-0.1148 


-0.1166 


-0.1155 


0.0021 




z 


0.2450 


0.2437 


0.2447 








0.2385 




0.2427 


0.2436 


0.21*1*4 


0.2365 


0.0042 




2n 


0.00287 
0.00149 












0.00027 


0.00200 


-0.00048 


-0.00061 




-0.00124 


0.00163 




§22 

B 33 
R 12 

P 23 












0.00238 




0.00261 






0.00162 
0.01944 


0.00180 




0.01905 












0.02399 




0.03084 


0.03061 




0.00487 




0.00066 












-0.00127 




-0.00077 


-0.00099 




-0.00050 


0.00113 




-0.00081 












0.00014 




0.00118 


0.00117 




0.00159 


0.00176 




-0.00057 












-0.00069 










-0.00363 


0.00216 


Atom 4 


X 


-0.4654 


-0.4682 


-0.4701 








-0.4702 








-0.4688 


-0.1*717 


0.0021 


(IX) 


y 


0.3272 


0.3259 


0.3258 








0.3255 








0.3263 


0.321*5 


0.0016 




z 


-0.0520 


-0.0573 


-0.0571 








-0.0629 




-0.0549 


-0.0552 


-0.0581 


-0.0515 


0.0039 




Si' 

fi 22 
pa 

« 12 

e 13 
B 23 


0.00266 












0.00357 




0.00522 


0.00522 




0.00403 


0.00166 




0.00120 












0.00070 




0.00379 


0.00388 




0.00323 


0.00146 




0.01613 












0.01697 




0.00943 


0.00926 




0. 00534 


0.00401 




0.00127 












0.00089 




0.00271 


0.00272 




0.00218 


0.00107 




-0.00122 












0.00017 




-0.00212 


-0.00213 




-0.00128 


0.00172 




0.00174 












0.00193 










0.00194 


0.00124 


Atom 5 


X 


0.4647 


0.4669 


O.U661 








0.4626 




0.1*639 


0.1*637 


0.1*634 


0.4630 


0.0022 


(X) 


y 


-0.3391 


-0.3377 


-0.3373 








-0.3369 




-0.3377 


-0.3377 


-0.3374 


-0.3372 


0.0017 




z 


0.175S 


0.1748 


0.1759 








0.1735 




0.1784 


0.1785 


0.1753 


0.1798 


0.0044 




S22 

a 33 

r' 2 

a 13 
P 23 


0.00200 












0.00383 




0.00285 


0.00278 




0.00194 


0.00187 




0.00142 












0.00269 




0.00048 


0.00039 




-0.00055 


0.00180 




0.01636 












0.01668 










0.01066 


0.00465 




0.00135 












0.00189 




0.00075 


0.00071 




0.00029 


0.00131 




0.00007 












-0.00176 




0.00105 


0.00097 




0.00029 


0.00190 




-0.00150 












0.00049 




-0.00002 


0.00000 




-0.00090 


0.00143 


Atom 6 


X 


0.1132 


0.1192 


0.1132 0.1218 


0.1130 0.1178 


0.1130 


0.1163 


0.1172 




0.1158 


0.1157 


0.1139 


0.1159 


0.0026 


N (1) 


z 


0.5065 


0.4959 


0.5000 0.5043 


0.5100 




0.5055 


0.5057 




0.5110 


0.5111 


0.5100 


0.5110 


0.0044 




a 33 

a 12 

D 13 


0.0038 






0.00791 




0.00634 


0.00600 




0.00658 






0.00488 


0.00295 




0.01480 






0.01356 




0.01021 


0.01008 




0.00970 






0.00866 


0.00430 




0.0032C 






0.00934 


0.00600 


"0.0 1036 


0.01189 


0.00500 


0.00792 


0.00711 




0.00646 


0.00412 




0.00000 






-0.01098 


-0.00500 


-0.00476 


-0. 00*153 










-0.00307 


0.00432 


Atom 7 


X 


-0.1304 


-0.1284 


-0.1284 








-0.1340 




-0.1306 


-0.1309 


-0.1321 


-0.1317 


0.0025 


(VIII) 


Z 


-0.1208 


-0.1209 


-0. 1212 








-0.1187 




-0.1155 


-0.1156 


-0.1165 


-0.1155 


0.0053 




n 33 

B 12 
I3 


0.00649 












0.00105 


0.00400 


0.00484 


0.00501 




0.00464 


0. 00265 




0.01200 












0.01251 










0.00971 


0.00512 
0.00354 




0.00514 












-0.00147 


0.00300 


0.00361 


0.00376 




0.00401 




-0.00700 












-0.00487 




-0.01211 


-0.01208 




-0.01273 


0.00330 


Atom 8 


X 


0.1426 


0.1366 


0.1333 








0.1333 








0.1347 


0.1330 


0.0023 


(VI 1) 


z 


0.1175 


0.1184 


0.1185 








0.1123 
0.00459 




0.1180 


0.1176 


0.1169 


0.1198 


0.0039 




n 33 

a 12 
P I3 


0.00728 














0.00087 


0.00071 




0.00014 


0.00228 




0.00710 












0.00457 




0.01211 


0.01 194 




0.01199 


0.00309 




0.00710 












0.00414 




0.00045 


0.00021 




-0.00011 


0.00304 




-0.00545 












-0.00608 




-0.00294 


-0.00276 




-0.00440 


0.00292 


Atom 9 


X 


-0.2912 


-0.2927 


-0.2908 








-0.2943 




-0.2941 




-0.2933 


-0.2935 


0.0020 


(IV) 


a 33 

a 12 


0.4756 
0.00823 


0.4738 


0.4762 








O.U723 
0.00792 




0.4755 


0.4755 


0.4725 


0.4814 
0.00342 


0.0039 
0.00311 




0.00784 












0.00919 




0.00748 


0.00731 




0.00627 


0.00324 




0.00627 












0.00788 




0.00688 


0.00686 




0.00192 


0.00398 




0.00209 












0.00155 










-0.00473 


0.00327 


Atom 10 


X 


0.3699 


0.3707 


0.3695 


0.3692 




0.3707 


0.3741 




0.3734 


0.3736 


0.3736 


0.3755 


0.0025 


(III) 


Z 


-0.4050 


-0.3958 


-0.3940 


-0.1*050 -0.3928 




-0.3946 


-0.397'* 




-0.3919 


-0.3921 


-0.3954 


-0.3878 


0.0043 




1" 

ft 33 

13 


0.00651 






0.00743 




0.00772 


0.01052 




0.01176 


0.1206 




0.00812 


0.00357 




0.01627 






0.01411 




0.01235 


0.01374 










0.00369 


0.00462 




0.00640 






0.00495 




0.00596 


0.00862 




0.01144 


0.1182 




0.00809 


0.00467 




0.00233 






0.00652 




0.00401 


0.00653 










-0.00179 


0.00437 


Atom 11 


X 


-O.'lU 12 


-0.4433 


-0.4428 








-0.4402 




-0.4444 


-0.4445 


-0. 4l» 75 


-0.1*1*50 


0.0023 
0.0041 


(II) 


z 


0.2764 


0.2773 


0.2752 








0.2821 




0.2822 




0.2821 


0.2803 




f,33 
•12 
> 5 I3 


0.00200 












-0.00019 


0.00200 


0.00024 


0.00009 




0.00067 


0.00180 




0.01777 












0.00976 




0.01433 


0.0129U 




0.01353 


0.00356 




-0.00009 












-0.00144 




0.00209 


0.00165 




0.00292 


0.00315 




0.00400 












-0.00003 




0.00245 


0.00216 




0.00258 


0.00319 


Atom 12 


X 


0.4538 


0.4553 


0.4533 








0.1*581 




0.1*567 


0.1*568 


0.U551 


0.1*566 


0.0021 


(1) 


s» 

a 33 

a' 2 
P 13 


-0.1611 

0.00200 


-0.1629 


-0.1625 








-0.1669 
0.00283 




-0.1632 
0.00427 


-0.1628 
0.00424 


-0.1607 


-0.1674 
0.00473 


0.0040 
0.00241 




0.00910 












0.01021 




0.00641 


0.00666 




0.00790 


0.00369 




0.00000 












0.00377 




0.00603 


0.00593 




0.00677 


0.00355 




0.00200 












0.00134 










0.00599 


0.00404 


Atom 13 


X 


0.3479 


0.3476 


0.3477 








0.3477 








0.3474 


0.3474 


0.0008 


S (1) 


z 

a" 
a 33 

a 12 

p 13 


-0.2433 
0.00163 


•0.2447 


-0.2428 








-0.2414 
0.00174 




-0.2433 


-0.2432 


-0.2423 


-0.2450 
0.00118 


0.0030 
0.00077 




0.00714 












0.00684 




0.00661 


0.00655 




0.00502 


0.00126 




0.00067 












0.00158 




0.00155 






0.00126 


0.00105 




-0.00041 












-0.00229 




-0.00141 


-0.00148 




0.00003 


0.00099 


Atom 14 


X 


-0.3174 


-0.3189 


-0.3187 








-0.3187 








-0.3183 


-0.3189 


0.0008 


S (II) 


«33 

a' 2 

P 13 


0.3144 
0.00194 


0.3123 


0.3137 








0.3156 
0.00223 




0.3135 
0.00215 


0.3136 


0.3147 


0.3112 
0.00172 


0.0030 
0.00087 




0.00756 












0.00871 




0.00982 


0.00987 




0.00558 


0.00181 




0.00082 












0.00209 




0.00179 


0.00186 




0.00198 


0.00115 




0. 00017 












-0.00142 




-0.00048 


-0.00048 




0.00069 


0.00106 


Atom 15 


z 


0.5350 


0.5529 


0.5467 








0.5338 


0.5500 


0.5535 


0.5525 


0.5432 


0.5520 


0.0091 


C (II) 


P 33 


0.00238 












0.00371 




0.00343 


0.00336 




0.00208 


0.00105 




0.01029 












0.00672 




0.00244 


0.00236 




-0.00152 


0.00364 


Atom 16 


z 


0.0572 


0.0571 


0.0575 








0.0573 








0.0574 


0.0587 


0.0018 


A* (II) 


P 33 


0.00189 












O.OC089 




0.00102 






0.00022 


0.00026 




0.00770 












0.01072 




0.00986 


0.01009 




0.00451 


0.00235 


Atom 17 


z 


0.5072 


0.5050 


0.5060 








0.4887 


0.5000 


0.5015 




0.1*948 


0.4994 


0.0088 


C (1) 


8 33 


0.01033 












0.01027 




0.02050 


0.02009 




0.01523 


0.00689 


Atom 18 


33 


0.00891 












0.00630 




C. 00614 






0.00511 


. 002 1 1 


A-t (1) 




























R 




0.20U 




0.176 


0.177 


0.167 






0.167 












Weighted 




0.149 




0.119 


0.117 


0.102 






0.097 












s 




2.19 




1.05 


1.90 


1.69 






1.60 












Est imate 


d s 




1.90 


1.86 


1.79 




1.68 


1.59 




1.50 


1.1*8 


1.58 


1.1*0 




No. of d 
in L.S. 
calculat 


lata 




534 


496 


563 




546 


546 




546 


546 


546 


546 




i on 




























Tempera t 


ure 








6 




6 


2,3,6 




3,6,7 
3,11,12 


3,6,7 




2,3,5 




factor not 














7,3,11 




8,11,12 




6,7,8 




pos i t i ve 
















12 










9,11,12 




def i n i te 


























13,14,15 
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Table V — Some Interatomic Distances Obtained from 
Least Squares Calculations (Second Set of Weights) 



Distance 


Cycle 1' 
A 


Cycle 2' 
A 


Ctf)-N(l) 


1.40 


1.43 


C(II)— N(II) 


1.29 


1.25 


S(I)-()(V) 


1.46 


1.44 


Sd)-O(III) 


1.38 


1.38 


Sd)-O(I) 


1.46 


1.44 


S(II)— O(VI) 


1.47 


1.48 


S(II)-OdV) 


1.48 


1.50 


8(H)— (II) 


1.50 


1.49 


A1(I)— O(VII) 


1.92 


1.89 


Al(I)— O(VIII) 


1.86 


1.86 


Al(II)— O(IX) 


1.90 


1.92 


Al(II)— 0(X) 


1.91 


1.91 



In cycle 2' the starting parameters were the same as those resulting 
from cycle 1' (new weights) except for the .r-parameters of N(I) and 
N(II) and the z-parameter of N(I). Also, it was found that under the 
conditions set for the weighting in cycle 1', only 496 amplitudes should 
have been weighted unity. For these reflections and the starting parame- 
ters shown in Table III, the R value was 0.176, weighted R, 0.119 and 
error of fit, 1 .85. Again only scale and positional parameters were allowed 
to vary. Changes were not large except for the N and C(II) parameters. 
Some distances calculated from these parameters are given in Table V. 
(C — N distances are always on assumption of planarity of the guani- 
dinium group.) Note that again the C(I) — N(I) distance is short, the 
C(H) -N(II) long, but the average is the expected value for such a 
bond. Also noteworthy is the continued tendency of S(I) — O(III) to 
be short. In fact, there is a tendency throughout for the S(I) — O dis- 
tances to be shorter on the average than the S(II) — O distances. Ex- 
amining the correlation matrix for this cycle we may summarize the 
results as follows (Table VI). Only those pairs for which | p | ^ 0.40 are 
listed. Thus of the 946 p u (i 5* j) terms only 75 are ^0.40. Important 
also is the fact that a large number, 677, of the terms are less than 0.10, 
many much less than 0.10; 194 of the | p (J | lie between 0.10 and 0.40. 
These could be important especially if one parameter has many inter- 
actions of moderate size with other parameters. 

Earlier we gave some examples of | p,y | that were calculated from the 
inverse matrix of cycle 12 (old weights). It is seen from examination of 
Table \T that the values for the particular | p - | obtained from cycle 2' 
arc essentially the same except for the value for the x N (ii) — %d) in- 
teraction. The value is much higher, 0.62, than the one, 0.10, obtained 
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from the incorrect orientation of the guanidinium ions. Thus, incorrect 
values for parameters can uncouple parameters. Furthermore, this ap- 
pears to be the reason that there was not much change in the incorrect 
//-parameter of X(II) in the first three cycles. That is, a parameter which 
is given a value which tends to make it independent may not change 
rapidly to a value which tends to make it dependent. 

The purpose of the next cycle was to see the results of allowing the 
parameters, both positional and thermal, of only the N and 0(111 ) 
atoms to vary. Before carrying out this calculation, however, the posi- 
tions of hydrogen atoms were estimated. The guanidinium ions were 
considered to be essentially planar, and the z-parameters of the guanidi- 
nium H's taken as 0.55 for those about the threefold axes at §,§ and 
§ ,$, and 0.505 for those about the axis at 0,0. For the water molecules, 
the links with the S() 4 oxygen atoms were considered and the tilt of the 
water molecule estimated accordingly. In any given level of H.O mole- 
cules about either of the nonequivalent axes, the z's were taken equal. 
The H— O— H angle was taken as 105° and the O— H distance, 0.9(5 
A. (The initial H-parameter values will not be listed; however, the 
last set used will be listed later.) First, H contributions to the F h k-i for 
h,k,l positive were calculated for two different orientations of the 
guanidinium ions, namely: ( + . — ) and ( + , + ). (The program used 
for this calculation was written by R. G. Treuting; the atomic scattering 
factors for H were those of Viervoll and Ogrim. 16 ) These calculations, 
together with consideration of previous calculations of the amplitudes, 
corroborated the conclusion that the (+,— ) orientation was the most 
probable one. 

The N-parameters were readjusted to yield the most probable C — N 
distance, and the z-parameter of 0(111 ) was started at —0.405. Those 
observed amplitudes with sin 2 0/A 2 < 0.0800 which were not strongly 
affected by extinction were rewcighted unity. The total number of re- 
flections weighted unity was 5(58. The H atoms were put into the calcu- 
lation as "fixed atoms" (see Ref. 13) with isotropic temperature factor 
B = 3.00 A". The over-all R value was 0.177, weighted R, 0.117, and 
error of fit, 1.90. 

The results of the least squares calculation are given in Table IV 
cycle 3'. It is seen that the O(III) z-parameter returned to that of the 
previous cycle. The X(I) .r-parameter increased somewhat, implying a 
C(I) — N(I) distance of 1.37 A. The parameters of X(II) imply a 
C(II)— N(II) distance of 1.33 A. 

In Table VII, we list those correlation coefficients greater than or 
equal to 0.40. If this table is compared with Table VI, one finds that the 
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coupling of N(I) and N(II) positional parameters is still as strong as in 
the previous cycle. In both cycles 2' and 3', the correlation matrices 
showed no strong interaction between O(III) and nitrogen atom 
parameters. The correlation matrix of cycle 3' indicated that there are 
some very strong interactions in pairs of thermal parameters. As ex- 
pected, there was corroboration of a strong interaction between the 
/3 33 's of the N atoms. 

For this case, it might be worthwhile to show the Ap.'s obtained from 
the complete solution of the 21 normal equations compared with those 
obtained from the diagonal term approximation. These are given in 
Table VIII together with the a's calculated by the Busing-Levy pro- 
gram. As expected, several of the Ap/s for particular i are quite different, 
particularly for those which are highly correlated (see Table VII). 

Before proceeding to the next cycle, the calculated and observed data 
were examined for any outstanding discrepancies and rechecks were 
made on the intensity data. It was found that 27 of the reflections which 
were listed as observed should have been listed as unobserved. It was 
also found that 5 reflections which were recorded as unobserved should 
have been observed by the instrument but were missed. These -were 
obtained from film data. 

Slight changes were made in the H-parameters; the ^-parameter of 
N(I) was returned to 0.113 and necessary changes made in the 19 and 
0i3 thermal parameters of N(I). Now the Busing-Levy program calcu- 
lates and stores all derivatives, so that it is possible to allow different 
sets of parameters to remain constant and solve for sets of Ap, for each 
initial set of parameters. In cycle 4'a, therefore, we first allowed only 
the N(I), N(II), and O(III) parameters to vary and then in 4'b, 



Table VII — Correlation Coefficients from Cycle 3' 



IpI 


tj 


0.40-0.50 


1,2 


7,12 


8,9 


8,15 


0.50-0. GO 


4,7 








0.60-0.70 


1,10 


5,7 


5,14 




0.70-0.80 


3,11 


0,13 






0.80-0.90 








■- 


0.90-1.00 


12,14 


18,20 







Parameter numbers 





X 


y 


i 


0ii 


fin 


033 


01* 


flu 


0» 


N(II) 


1 


2 


3 


4 


5 


<) 


7 


8 


9 


N(I) 


10 




11 


12 




13 


14 


15 




O(III) 


10 




17 


18 




19 


20 


21 
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Table VIII — Parameter Changes and Error Estimates 
from Cycle 3' 



Parameter number 


Busing-Levy 


Diagonal term 


o-'s from 


(see Table VII) 


approx. 


Busing-Levy 


1 


-0.0029 


-0.0020 


0.0027 


2 


-0.0050 


-0.0027 


0.0022 


3 


0.0058 


0.0046 


0.0037 


4 


-0. 00601 


-0.00105 


0.00268 


5 


0.00270 


0.00418 


0.00245 





0.00231 


0.00477 


0.00420 


7 


-0.00273 


-0.00136 


0.00202 


8 


-0.00227 


-0.00333 


0.00271 


9 


-0.00371 


-0.00053 


0.00225 


10 


0.0048 


0.0051 


0.0025 


11 


0.0057 


-0.0030 


0.0048 


12 


0.00411 


0.00023 


0.00290 


13 


-0.00124 


0.00172 


0.00527 


14 


0.00014 


0.00113 


0.00400 


15 


-0.01098 


-0.00983 


0.00552 


16 


-0.0003 


-0.0009 


0.0014 


17 


0.0121 


0.0120 


0.0023 


18 


0.00092 


0.00160 


0.00228 


19 


-0.00210 


-0.00133 


0.00311 


20 


-0.00145 


-0.00230 


0.00307 


21 


0.00419 


0.00025 


0.00351 



varied all parameters except the scale factors. The results are shown in 
Table IV. Again in both cases, the N(I) .r-parameter increased; there 
were changes in the N(II) parameters, but the implied C — N distance 
1.35 A was good. Also the ^-parameter of O(III) seemed to improve, 
especially when all the parameters were allowed to vary. But in 4'a, 
the thermal parameter matrix of the N(I) atom was not positive defi- 
nite, while in 4'b, seven atoms had thermal parameter matrices which 
were not positive definite. Also there were continued oscillations and 
large error estimates. It was evident that real convergence would not be 
attained. 

However, because the N and 0(111 ) parameters did look encouraging, 
it was decided to try one more cycle. This time the parameters of the 
water hydrogen atoms were recalculated in a somewhat different way. 
In a recent paper, IJ Aleksandrov, Lundin and Mikhailov report results 
of a study of the distribution of hydrogen atoms in guanidinium alumi- 
num sulfate hcxahydrate by means of proton magnetic resonance experi- 
ments. They report that the nearest neighbor p — p (proton-proton) 
vectors are perpendicular to the ct\ , a 2 and a 3 axes.f They argue that on 
the basis of symmetry considerations all II atoms bonded to O's in a 

t Previously, Spence and Muller 18 had reported this to be so for the p — p 
vectors of the water molecules, but had concluded that the p — p vectors of the 
guanidinium groups could be parallel to the c-axis with a separation of 2.05 A. 
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single octahedron layer about a threefold axis must have the same 
2-parameter. Of course, this is true only for those hydrogen atoms 
bonded to N(I) atoms and to the water molecules about the threefold 
axis at 0,0. The trigonal axes and planes of symmetry are such that only 
three atoms about the axis at |,f and three about the axis at § ,f must 
have the same value of z. 

Thus contrary to the statements of Aleksandrov et erf, 19 symmetry 
conditions do not require all the nearest neighbor H — H vectors to be 
parallel to the (00-1) plane, nor must they all be perpendicular to the 
ai , a 2 and a :) axes. Only for those about the threefold axis where the 
mirror planes intersect, namely at 0,0 must this be the case. However, 
it is possible that the nearest neighbor H — H vectors about the three- 
fold axes at 3 f, f f are close to parallelism with the (00-1) plane and 
perpendicularity to the ai , a 2 , a A axes. 

Furthermore, Aleksandrov et al 19 refer to the trial structure reported 
by Varfolomeeva et al. 3 Although that structure is incorrect, it would 
have no noticeable effect on the conclusions of Aleksandrov et al, since 
they discuss only the nearest neighbor H — H vectors. 

Thus, in calculating the H parameters, the tilting of the water H — H 
bonds out of the (00- 1 ) plane and skewness to the a\ , a 2 , a 3 axes was 
permitted in those water molecules about the threefold axes at $ $, f $. 
(The guanidinium ions, however, were assumed to be planar.) In calcu- 
lating the H positions, the water molecules were assumed to lie in the 
planes connecting the water oxygen atom with the two sulfate oxygen 
atoms involved in the hydrogen bonding. The bisector of the H — O — H 
angle of 105° was taken as the line passing through the center of the 
water oxygen atom and the center of the line connecting the two sul- 
fate oxygen atoms involved. The parameters of the N and O atoms 
involved were those from cycle 4'b. The H-parameters thus deduced are 
listed in Table IX. The new parameters caused some differences in the 



Table IX — H Parameters Used in Final Cycle 



Description 


X 


y 


% 


on N (I) (atom 6) 


0.205 


0.086 


0.51 


on N (II) (atom 1) 


0.465 


0.256 


0.56 




0.564 


0.434 


0.56 


on O (VIII) (atom 7) 


0.139 


0.218 


-0.148 


on 0( VI I) (a torn 8) 


-0.072 


0.134 


0.156 


on 0( IX) (atom 4) 


0.457 


0.257 


-0.124 




0.526 


0.400 


-0.111 


on O(X) (atom 5) 


-0.452 


-0.260 


0.205 




0.464 


0.588 


0.219 
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contributions to several amplitudes, but in general not very important 
ones. 

Some necessary adjustments of thermal parameters resulting from 
cycle 4'b were made. In cycle 5'a,t only those positional parameters 
were varied in which changes greater than a/o occurred between previ- 
ous cycles 2' and 4', all thermal parameters were varied in which there 
were changes greater than a/5 between cycles 1' and 4'; all scale factors 
were kept constant. In 5'b, only those parameters were varied in which 
changes in 5'a were greater than a/5. In 5'c, only positional parameters 
were varied. In 5'd all parameters were varied. All results are listed in 
Table IV. Differences range from very small to very large and are in- 
dicative of the unattainability of convergence. We list also the a'sl in 
the Ap.'s for the last cycle 5'd in the last column of Table IV. These 
are especially large for most of the thermal parameters and for most of 
the z-parameters, and reflect the strong interdependence in pairs of 
parameters. 

The correlation matrixt for cycle 5'd contains 6,670 p#(* ^ j) terms. 
Thus we shall again only list the values of | pa | ^ 0.40 (Table X). Of 
the 6,670 terms in the matrix, 176 have values greater than 0.40; 1,889 
have values greater than 0.10. 

On examining Table X, one finds that no interactions of scale factors 
with positional parameters are listed. In fact, the correlation coefficients 
for such combinations are all very low. However, there are all the other 
types of interactions, namely: scale factor-thermal parameter, thermal 
parameter- thermal parameter, positional parameter-positional parame- 
ter, and several (those with asterisk) positional parameter-thermal 
parameter. Most often, also, the interdependence is between analogous 
parameters; for example, a ^-parameter of an atom interacts with z- 
parameters of other atoms. Even when a positional parameter inter- 
acts with a thermal parameter, an analogy exists, e.g., a ^-parameter 
interacts with a fo-parameter. This makes physical sense, of course, 
and gives us some confidence that the correlation coefficients reflect the 
structural interdependence of the parameters. Correlation may be 
caused partially by the experimental teehnique§ but it is unlikely to 
result mainly from the ill-conditioning of the normal equations by a 

t It, Should be kept in mind that all cycles 5' refer to the derivatives evaluated 
with the parameters of cycle 5'a. 

% It is worth emphasizing that statistical theory precludes the use of the error 
estimates or normal ecpiations matrix to determine the statistical significance of 
the parameters listed. Only if convergence is actually attained can these numbers 
be so used. Nevertheless, in a practical way, the error estimates and correlation co- 
efficients do give us important information in the course of refinement or, as in the 
present case, relative to the unattainability of convergence. 

§ X-ray vs neutron diffraction. 
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Table X — Correlation Coefficients From Cycle 5'd.f 
(Only I pa \ > 0.40 Are Listed) 



* — jitji.jii 



0.40-0.50 



0.50-0. GO 

0.G0-0.70 

0.70-0.80 

0.80-0.90 
0.90-1.00 



3-4,5; 4-105; 6-75; 7-75; 8-75; 11-38,47,61,67,85; 15-56; 18-19,34 
20-47 J,61, 97, 103; 21-30; 24-30,86,88; 25-26,27; 29-47,61; 30-33 
36-37; 38-61,103; 40-51,62; 42-68; 43-66$; 44-46$; 45-46; 47-55 
67,73,85,91,103; 51-64; 52-53; 55-97,103; 57-59; 61-73,79,85,91 
63-69; 67-103; 75-81 ; 77-79$,83; 79-83$,97,103; 81-83; 97-108,114 
103-108,114; 105-113; 108-110$; 113-116 

4-5,6,7,8; 5-113; 8-99,113; 11-97,103; 12-15; 13-15,58; 16-59; 
19-84; 20-91; 22-24; 26-89; 27-941; 29-85,97,103; 30-94; 37-53$; 
38-47,97; 39-42; 41-50,63,69; 43-71; 47-97; 48-51; 50-69; 51-62; 
61-97,103; 67-97; 73-83f,97,103; 85-91,97,103; 110-115 

5-105; 13-56; 14-57; 20-29; 21-24,88; 23-86; 27-92$; 28-90; 30-92; 
37-46; 40-42,49,68,70; 49-64; 52-65; 72-78,801,82$; 74-78$; 76- 
78$; 86-88; 96-104$, 106$; 98-102$; 100-102$; 108-115$; 111-116$ 

5-6,7,8; 8-105; 9-54; 18-27; 21-86; 32-93; 35-95; 36-66; 45-60; 
49-51,62; 50-63; 96-102; 108-114; 111-113$ 

6-7,8; 11-55; 25-91; 38-67; 62-64; 73-79; 98-100 

7-8; 47-61; 56-58; 68-70; 74-70; 80-82; 92-94; 97-103; 104-106 



t See last column of Table IV for parameter numbers. 
$ Positional -thermal parameter correlation. 

reasonable but not necessarily ideal weighting technique. It will be 
noticed also that the same pairs of parameters show very nearly the 
same measure of interdependence as indicated by earlier calculations, 
again corroborating the point that it is the structural model (including 
atomic form factors) which causes the interactions. 

For the sake of completeness, we show in Table XI a list of observed 
amplitudes compared with those calculated from the parameters used 
initially in cycle 5' and including the contributions of the H atoms with 
parameters shown in Table IX. Including consideration of multiplicity 
and the differences when calculated amplitudes are greater than the 
threshold values (with half the threshold value included in the denom- 
inator) for reflections not observed, the discrepancy factor is 0.11. t 

The over-all agreement is quite good despite several discrepancies in 
which a calculated amplitude is above the threshold value for an unob- 
served reflection.! Table XI attests to the validity of the conclusion 
that the general features of the structure are correct. 

t Six amplitudes, those of reflections 30-0, 111, 21-1, 221, 42-1 and 21-2, suffer- 
ing from extinction were excluded in calculation of this discrepancy factor. 

$ These are a product of the instrument which sometimes missed reflections, 
which, according to visual estimates of photographic intensities, it should not have 
missed. 
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VI. FURTHER COMMENTS ON THE INDETERMINACY OF THE EXACT STRUC- 
TURE OF GUANIDINIUM ALUMINUM SULFATE HEXAHYDRATE 

0. l Importance of the Weighting Procedure 

The use of two very different weighting procedures did not break down 
the high correlations existing between parameters. It is doubtful, 
especially in the case of so large a number of parameters, that any 
reasonable weighting procedure would succeed in uncoupling the parame- 
ters sufficiently to lead to greater determinacy. 

fi.2 Effect of Keeping Some of the Parameters Constant while Alloiving 
Others to Vary 

In the case that there is correlation between parameters, it would seem 
that, at least in the final stages of the refinement, holding of such parame- 
ters constant could lead to erroneous results. In a case involving a smaller 
number of parameters it might be possible to obtain a confidence region 5 
for all the parameters by holding some of the parameters constant, but 
at several different values. For example, suppose the problem involves n 
almost independent parameters and two almost completely dependent 
parameters which appear to prevent convergence. Choosing several 
judicious values of one of the latter and making the calculation for each 
one will give sets of values for the other parameters which will allow the 
construction of the equiprobability ellipsoids. 

However, in a problem involving many parameters, and many large 
and multiple correlations, such a technique would appear to be im- 
practical. It should be mentioned that if the model were very nearly 
linear, only those correlations very near ±1 would be important in the 
unattainability of convergence. However, it is possible that the more 
nonlinear the model, the more important the other correlations become. 

0.3 Possible Effects of Increasing the Number of Observed Data 

There are two ways in which the number of data might be increased. 
One is to obtain more of the weak intensities by increasing the detector 
sensitivity. It does not seem that this would have the effect of decreas- 
ing the correlations. This was shown to some extent by the calculations 
based on the two different weighting schemes. In the first case the 
weighted evaluated derivatives for unobserved reflections were included; 
in the second, these were given zero weight and therefore excluded. 
Also, the exclusion of reflections for which sin 2 0/X 2 < 0.0800 did not 
have an apparently significant effect on the correlations. (Compare, for 
example, analogous values in Tables VI and X.) 
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The other way in which to increase the number of data is to use shorter 
wavelength radiation. Now, it is not necessary actually to measure these 
data before determining the effect on the correlations because the cor- 
relation coefficients, as calculated, depend only on the model and the 
evaluated derivatives. It is unlikely that the situation would change very 
much if the additional terms were included because the relationship of 
the derivatives with respect to correlated parameters would probably 
not change very much. 

In the case of tetragonal BaTi0 3 , M higher index reflections would 
have almost no important contributions from the oxygen atoms. Thus 
the interactions among oxygen ion parameters will not be affected. 
Similarly, interactions among the metal ion parameters will probably 
not be much affected. But interactions between the two groups could be 
reduced. However, in the case of an all light atom structure, it would 
appear that the extra data would probably not reduce the correlations. 

6.4 Possible Effect of Greater Accuracy in Measurement of Observed In- 
tensities 

The effect of greater accuracy in measurement of the observed in- 
tensities is not really predictable in this case. To be sure, in each iteration 
the reduction of s = V 2(\/wA)' 2 /\/ni — n would reduce the apparent 
size of the equiprobability surfaces. This we certainly know. 

However, we must ask first whether there is a limit to the accuracy of 
the observed amplitude. One would suspect that there is such a limit. 
Furthermore as pointed out by Caticha-Ellis and Rimsky, 2 there will 
always be a discrepancy between the calculated and true values of the 
amplitudes. Thus s has a lower positive limit. 

Reduction of s would not only decrease the size of the equiprobability 
surfaces (and therefore, of course, the standard estimates of error) Jmt it 
would also decrease the components of the vector v, Vj = 2(\/W)j)- 
(\/w>A), where the Z), are the evaluated derivatives. Thus, for example, 
if cycle o'd were repeated with each A decreased to £ of its value, each v, , 
and therefore each Ap, = ^2j &,-/«,- would be reduced to the same extent. 
Of course an average reduction of \ might not do the same thing. In fact, 
with a poor distribution of the reduction in A, the Ap, in some cases 
could even be larger, depending on the algebraic values of the D, . 

Actually the nature of the shape of the equiprobability surfaces might 
give the best clue to what might happen if increased accuracy of measure- 
ment were attainable. The nonlinearity of the model would probably 
play an important part. The more nonlinear, the more important are apt 
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to be those correlations which are not perfect. Of course, even one perfect 
correlation ±1, renders the whole problem indeterminate 5 if insistence 
is made on allowing all parameters to vary in an iteration. This is not 
necessary, however, and one could learn a great deal about the parame- 
ters of a structure which has only one perfect correlation and the rest 
very small ones (see Section 6.2). In the present case, there are many 
correlations having absolute values between 0.90 and 1.00 (Table X). 
These have the specific values: 0.917, 0.905, 0.913, 0.907, 0.975, 0.963, 
0.901, 0.979, and 0.902, respectively. Perhaps the most important ones 
are the three closest to unity. 

In the case of gross nonlinearity it seems possible that these and so 
many of the other high correlations of Table X could cause unattain- 
ability of convergence even if the lowest limit of s were attained. That is, 
the shape of the equiprobability surface may be such as to prevent the 
practical attainment of separate estimates of the parameters (see also 
Ref. 21) from the given data. This seems to be true of the BaTi0 3 

5,6 

case. 

Needless to say, a measure of doubt remains. Further work might aid 
in removing this doubt. This would involve trying to obtain more data 
and of greater accuracy, and further calculations. Our doing this is not 
presently contemplated. 

G.5 Fourier Synthesis vs Least Squares 

In the case of tetragonal barium titanate, Fourier synthesis produced 
no improvement on the least squares method." It is likely that with the 
present data, the situation in the case of the G.A.S.H. would be the 
same. On the other hand, there is no requirement of linearity in the 
Fourier synthesis: the actual amplitudes are the Fourier coefficients. 
In the least squares technique, an approximation is used: i.e., 



Fhki(pi,p< , • -,Pn) " Fiaiifh + &Pi,]h + Ap 2 , • • -,[)„ + Ap„; 

A/)j 



= F ll u(Fn,p„---,P, l ) + Z d ^ 

;-l dpj 



where [h ,fh , ■••,[)„ are approximate but nearly true values of the 
parameters. It is possible that higher order terms could be important 
here, but it is not clear that inclusion of the next higher order terms would 
necessarily lead to improvement. Also, the calculation would increase in 
complexity. 

Cochran has shown that a rather close relationship exists between the 
Fourier synthesis and least squares techniques. There are conditions on 
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this relationship given by Cochran 24 and Hoard and Geller 24 , and in 
addition in the actual least squares calculation, an approximation is 
made and nearness to linearity is assumed. Therefore, if the nonlinearity 
is not serious, convergence should be attainable in either case. If it is 
serious, the relationship could break down further and the Fourier 
.synthesis could conceivably converge when the least squares calculation 
tends not to converge. 

VII. COMMENTS ON THE SINGLE-CRYSTAL AUTOMATIC DIFFRACTOMETER 

As mentioned earlier, the data used in this work were collected four 
years ago. Since that time only one or two attempts were made to use 
the instrument for other studies. These were unsuccessful because of 
difficulties which are probably surmountable, but require modification 
of the instrument. 

The present instrument puts a lower limit on the sample size. To keep 
the time for recording a layer within reasonable bounds and to prevent 
the instrument from reacting to background scattering, only intensities 
above a certain preset count energize the circuitry which sets the crystal 
back and shifts speed. To obtain satisfactory counting rates the use of 
large crystals is required. (The intensity is proportional to the number 
of unit cells irradiated.) However, to obtain adequate or meaningful 
intensities from highly absorbing materials one must have small crystals. 
In short, the instrument presently is suited mainly to crystals with low 
absorption and from which sizable cylindrical specimens can be made. 

The indexing of the reflections was a tedious process. The possibility 
of error, particularly at the high angles, was great, but the use of photo- 
graphs and cross examination of data helped prevent errors. An improve- 
ment on the Bond-Benedict automatic single-crystal diffractometer 
would be provision for foolproof pre-indexing of the reflections. 

VIII. SUMMARY 

Extensive application of the least squares refinement technique 
(through the use of the Busing-Levy IBM 704 program) to three- 
dimensional X-ray data from crystals of guanidinium aluminum sulfate 
hexahydrate indicated that although the structure as originally reported 
for the isostructural guanidinium gallium sulfate is essentially correct, 
an exact structure is unattainable from the present data by means of the 
least squares method of refinement. The numerous high correlations of 
pairs of parameters, apparently linked with the nature of the structure, 
appear to be a primary cause of prevention of convergence. 
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The course of the calculations has been outlined with special emphasis 
on some of the more obvious parameter interactions, but tables are given 
to enable the more interested reader to examine the results in somewhat 
greater detail. 

The work also further demonstrates the importance of the correlation 
matrix as a tool for establishing the existence or nonexistence of inter- 
dependence of structural parameters. 
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